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Abstract 

Proton  exchange  membrane  fuel  cells  operating  with  Nation®  membranes  have  encountered  numerous  problems  associated  with  water 
management  and  CO  poisoning  because  of  their  low  temperature  of  operation.  Higher  temperature  membranes  have  been  investigated,  one 
such  membrane  being  polybenzimidazole  (PBI).  This  paper  presents  a  parametric  model,  which  predicts  the  polarization  performance  of  an 
intermediate  temperature  proton  exchange  membrane  fuel  cell  (PEMFC).  It  also  investigates  the  effects  of  porous  media  characteristics  on 
fuel  cell  performance.  Results  show  that  for  intermediate  temperature  fuel  cells,  mass  limitation  effects  are  absent  as  long  as  the  catalyst 
regions  are  sufficiently  permeable.  It  is  predicted  that  the  greatest  scope  for  improving  PBI  PEMFC  performance  is  increasing  the  membrane 
conductivity  and  improving  the  catalyst  performance,  as  it  interfaces  with  the  PBI  membrane. 

©  2005  Published  by  Elsevier  B.V. 
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1.  Introduction 

The  proton  exchange  membrane  fuel  cell  (PEMFC)  has 
long  been  considered  a  prime  candidate  for  future  energy 
systems.  However,  a  number  of  debilitating  factors — low 
electrode  kinetics,  catalyst  poisoning,  and  membrane  mate¬ 
rials  limitations — have  hindered  its  development. 

The  Nafion®  membrane  typically  used  in  low-temperature 
PEMFCs  requires  critical  water  management  for  sustained 
operation.  As  a  result,  the  fuel  cell  is  operated  below  90  °C. 
Proper  water  management  is  necessary  to  simultaneously  pre¬ 
vent  dehydration  of  the  membrane  and  flooding  of  the  elec¬ 
trodes.  This  critical  water  management  requirement  makes 
Nafion®  unsuitable  for  applications  with  variable  loading 
conditions,  such  as  automotive  applications. 

At  low  temperatures,  carbon  monoxide  (CO)  poisoning 
of  the  platinum  (Pt)  catalyst  is  also  a  problem.  Traces  of 
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CO,  present  in  the  hydrogen  feed,  are  preferentially  adsorbed 
on  catalyst  sites  resulting  in  reduced  efficiency  and  possible 
permanent  poisoning  of  the  catalyst. 

One  solution  to  these  problems  is  operating  the  PEMFC 
at  intermediate  temperatures  (150-200  °C).  At  such  tempera¬ 
tures,  CO  poisoning  becomes  less  prominent,  electrode  kinet¬ 
ics  are  faster,  and  water  would  exist  primarily  in  the  vapor 
phase  precluding  problems  associated  with  water  manage¬ 
ment  and  mass  transport  limitations.  However,  since  Nafion® 
requires  a  high  water  content  to  maintain  proton  conductivity, 
it  is  very  ineffective  at  temperatures  above  the  boiling  point 
of  water.  Therefore,  alternative  membranes,  which  maintain 
mechanical  strength  and  chemical  stability  at  elevated  tem¬ 
peratures,  are  needed. 

One  such  membrane  is  polybenzimidazole  (PBI).  PBI 
has  the  advantage  of  being  lightweight  as  well  as  stable 
at  higher  temperatures  (150-200  °C),  and  it  attains  proton 
conductivity  when  doped  with  phosphoric  or  sulfuric  acid 
[1-12].  Protons  are  transported  through  the  solid  matrix; 
so,  its  conductivity  is  less  dependent  on  water  content  than 
Nafion®.  Proton  conduction  in  PBI  may  be  based  on  the 
Grotthuss  (proton  hopping)  mechanism,  as  suggested  by 
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Nomenclature 

a  effective  surface  area 

c  specific  heat  capacity 

d  average  diameter  of  spherical  particles 

Dm  effective  gas  diffusivity 

/>[  effective  thermal  diffusivity 

E  potential 

/  weight  fraction  of  Pt  in  Pt/C 

F  Faraday  constant  (96,485  C  mol- 1 ) 

i  current  density 

if)  exchange  current  density 

j  reaction  rate 

k  thermal  conductivity 

Ap  permeability 

l  thickness 

m  mass  fraction,  mass  loading 

M  molar  mass 

P  pressure 

R  universal  gas  constant  (8.3143  J  mol-1  K-1) 

S  source,  entropy 

T  temperature 

u  mass-averaged  Darcy  velocity 

Greek  letters 

a  thermal  diffusivity,  charge  transfer  coefficient 

e  porosity 

y  concentration  parameter 

tp  electrical  or  ionic  potential 

v  kinematic  viscosity 

p  mass  density 

a  electrical  or  ionic  conductivity 


Bouchet  et  al.  [2,3],  PBI  membrane  thicknesses  are  lower  than 
Nafion®  because  of  its  higher  mechanical  strength.  A  reduc¬ 
tion  in  membrane  thickness  leads  to  lower  ionic  resistance, 
since  resistance  is  proportional  to  thickness.  In  addition,  PBI 
has  a  lower  permeability  to  hydrogen  and  methanol  than 
Nafion®,  resulting  in  less  crossover  of  fuel  at  elevated  tem¬ 
peratures. 

PBI’s  conductivity  can  be  further  enhanced  by  alloying 
with  inorganic  composites.  Several  studies  have  been  con¬ 
ducted  on  the  feasibility  of  using  PBI  and  its  composites  as 
membranes  in  PEMFCs  (Ref.  [4]  and  references  therein). 
Composites  of  PBI  have  been  developed  with  sulfonated 
polysulfones  (SPSF)  [5];  zirconium  phosphate  (ZrP),  phos- 
photungstic  acid  (PWA),  silicotungstic  acid  (SiWA)  [6];  inor¬ 
ganic  phosphomolybdic  acid  (PM012)  [7];  imidazole  and 
1 -methyl  imidazole  (Me-Im)  [8];  silicotungstic  acid  and  sil¬ 
ica  [9]. 

The  conductivity  of  PBI  and  its  composites  reportedly 
depend  on  the  temperature,  humidity,  doping  level  of  acid 
treatment,  and  method  of  preparation  of  the  membrane.  The 
conductivity  of  PBI  treated  with  phosphoric  acid  ranges  from 


Table  1 

Comparison  of  PBI  and  Nafion® 


tr(Sm  *) 

lm  (p-m) 

Total  resistance  (£2  cm2) 

PBI 

1-9 

40-80 

0.044-0.800 

Nafion® 

17 

230 

0.135 

0.2  to  6.8  S  m-1  [6,9-12].  Doped  with  sulfuric  acid,  the  con¬ 
ductivity  ranges  from  5  to  6  S  m-1  [13,14].  KOH-doped  PBI 
has  been  reported  to  have  an  ionic  conductivity  of  9  S  m- 1 
[13].  Among  the  composites  of  PBI,  one  with  0.15  wt.%  ZrP 
has  exhibited  a  conductivity  of  9.6  Sm-1  [6].  The  conduc¬ 
tivity  of  Nafion®  for  fully  humidified  membrane  conditions 
is  17  S  m-1  [15,16].  So,  it  is  evident  that  PBI-based  mem¬ 
branes  have  relatively  lower  conductivities  than  Nafion®. 
However,  since  PBI’s  permissible  membrane  thickness  is  less 
than  Nafion®,  its  ionic  resistance  is  of  the  same  order  of  mag¬ 
nitude.  Table  1  compares  typical  values  of  conductivity  and 
membrane  thickness  for  PBI-based  membranes  and  Nafion®. 
For  a  given  cross-sectional  area,  the  ionic  resistances  of  PBI 
and  Nafion®  are  in  the  same  order  of  magnitude. 

Much  work  has  been  done  in  PEMFC  modeling,  over 
the  past  15  years,  on  PEMFCs  using  Nafion®  membranes, 
i.e.  low-temperature  PEMFCs  [17],  Although  intermediate 
temperature  membranes  have  been  investigated  for  use  in 
PEMFCs,  to  date,  no  mathematical  model  has  been  pub¬ 
lished  for  any  intermediate  temperature  PEMFC.  This  paper 
presents  a  mathematical  model  of  a  PEMFC  working  with  a 
PBI  membrane,  or  a  PBI  fuel  cell  (PBIFC).  The  model  is  one¬ 
dimensional,  steady  state  and  non-isothermal.  It  attempts  to 
simulate  the  experimental  data  presented  by  Wang  et  al.  [18], 
and  also  to  determine  the  optimum  performance  of  a  PBIFC. 
Most  papers,  presenting  parametric  models,  focus  on  operat¬ 
ing  conditions.  This  paper  focuses  on  geometric  factors  such 
as  porous  media  characteristics,  and  membrane  and  catalyst 
properties. 


2.  Model  development 

2.1.  Assumptions 

For  intermediate  temperature  fuel  cells,  water  is  expected 
to  exist  only  in  the  vapor  phase,  therefore,  single-phase 
assumptions  apply.  The  gas  mixtures  are  assumed  to  behave 
ideally.  The  porous  media  and  catalyst  regions  are  assumed 
to  be  isotropic  and  macro-homogeneous,  and  the  membrane 
is  assumed  to  be  impermeable  to  gas  flow.  It  is  assumed 
that  gases  are  supplied  to  the  fuel  cell  at  a  sufficiently  high 
rate,  so  that  there  is  little  variation  in  fluid  properties  along 
the  gas  channels.  We  also  neglect  rib  effects.  Using  scaling 
arguments,  we  expect  the  most  significant  variations  in  trans¬ 
port  properties  to  occur  in  the  direction  perpendicular  to  the 
membrane  electrode  assembly  (MEA)  cross-section.  Thus,  a 
one-dimensional  treatment  of  the  MEA  is  adopted. 
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2.2.  Governing  equations 


The  following  conservation  equations  govern  the  phenom¬ 
ena.  The  Darcy  velocity,  u ,  the  source  term,  S,  the  tempera¬ 
ture,  P,  and  the  current  densities,  i  are  volume-averaged  over 
the  entire  porous  media  (solid  and  fluid  regions).  The  pres¬ 
sure,  P,  and  density,  p  are  volume-averaged  over  the  void 
region.  The  phase  potentials,  ip,  are  volume-averaged  over 
their  respective  phases  (solid  and  electrolyte  phases). 


V  ■  (pu)  =  Sm 

(i) 

sm  =  y  '  Si 

(2) 

The  continuity  Eq.  ( 1 )  has  a  non-zero  term  on  the  RHS ,  which 
represents  solid-fluid  phase  change  in  the  catalyst  layers.  For 
a  PBI  membrane,  H+  ions  are  transported  across  the  mem¬ 
brane  in  a  solid  state.  So,  there  is  a  loss  of  fluid  mass  at  the 
anode  catalyst  layer  and  a  gain  in  fluid  mass  at  the  cathode 
catalyst  layer.  Total  mass  is  always  conserved  at  every  point, 
and  overall  fluid  mass  in  conserved  since  the  H+,  formed  at 
the  anode,  re-converts  to  fluid  mass  at  the  cathode.  Across 
the  membrane,  the  velocity  is  zero. 

For  porous  media  in  fuel  cells,  the  momentum  equation 
scales  to  Darcy’s  law.  All  other  terms  (inertia,  body  and  vis¬ 
cous  terms)  are  negligible. 


VP  =  —  — —  u 
kp 


(3) 


Considering  that  a  non-conservative  form  of  the  continuity 
equation  is  used,  the  species  equations  must  be  re-derived. 
Conservation  of  individual  species,  i,  is  given  by. 


V  ■  (pju)  =  V  •  {pm ill)  =  V  ■  (pe DmV nij)  +  S,  (4) 


Expanding  the  LHS  and  combing  with  the  continuity  equa¬ 
tion. 


w,V  ■  (pu)  +  pu  ■  V(m,)  =  V  •  (peDmV/77;)  +  .S',  (5) 

pu  ■  V;n;  =  V  •  (peZ)mV/77;)  +  Si  —  miSm  (6) 

The  last  term  on  the  RHS  does  not  normally  appear  in  the 
species  equation  (actually,  it  is  normally  zero).  However,  for 
this  case,  where  a  non-conservative  continuity  equation  is 
used,  it  must  be  included.  If  Eq.  (6)  is  summed  for  all  species 
and  combined  with  Eq.  (2),  both  sides  of  the  equation  become 
zero. 

The  energy  equation  accounts  for  convection,  conduction 
and  heat  generation  due  to  ohmic  or  Joule  heating,  as  well  as 
heat  of  reaction. 


Catalyst  layer 


-X- 


Gas  Diffuser 


Membrane  proton 
conducting  region 


Active  catalyst  site 


\  Void  region 

Solid  electrically 


conducting  region 
Fig.  1.  Schematic  of  a  catalyst  site. 


v  •  »e  =  v  ■  (-crfvtpe)  =  j 


(10) 


/  „  \  Kreactant 

•ref/  /'reactant  \ 

J  —  a,0 


V  Pref  ) 


exp 


Pref 

aF 

PP 


exp 


aF 

- (0s 

PP 


0e) 


(0s  0e) 


(ID 


The  Butler- Volmer  Eq.  (11)  states  that  the  rate  of  charge 
transfer  is  driven  by  the  difference  in  potential  between  the 
two  respective  phases,  and  is  affected  by  the  concentration  of 
reactants  at  the  catalyst  sites.  The  reaction  rate,  /,  defines  the 
rate  of  transfer  of  solid-state  current  to  electrolyte  phase  cur¬ 
rent  in  the  anode  catalyst  layer,  and  vice  versa  in  the  cathode 
catalyst  layer.  Fig.  1  shows  a  schematic  of  the  catalyst  site 
illustrating  the  charge  transfers.  This  concept  of  two  phase 
potential  and  current  flow,  often  overlooked  by  most  fuel  cell 
modelers,  is  critical  to  fully  analyze  the  fuel  cell. 

The  exchange  current  density,  ip,  is  defined  in  terms  of  the 
active  catalyst  surface  area.  The  effective  surface  area,  a,  is 
defined  as  the  ratio  of  the  total  active  catalyst  surface  area  to 
the  total  catalyst  region  volume.  This  way  the  reaction  rate  is 
volume-averaged  over  the  entire  solid  and  void  phases. 


2.3.  Source  terms 

The  source  terms  are  given  in  terms  of  the  reaction  rates 
in  each  catalyst  layer. 


pcu  ■  VP  =  pcV  ■  (P>jVP)  +  St  (7) 

St  =  S0hm  +  Srxn  (8) 

The  charge-conservation  equations  in  the  electrolyte  and 
solid  phases  are  given  by, 

v  •  is  =  V  •  (-<ffV0s)  =  —  j  (9) 


So2  =  j 


Mo2 


AF 


Sh2o  =  —  j 


,4Th2o 


2  F 


Sh2  —  j 


•  mh2 


2  F 


(12) 

(13) 


(14) 
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5rxn  —  j 


^0e  0s 


(15) 


diffusivity  is  assumed  constant  at  a  given  temperature  and 
pressure,  and  independent  of  composition. 


•2  -2 

Is  le 

^"ohm  "P 


Teff 


(16)  2.5.  Boundary  conditions 


There  is  no  electrolyte  current  in  the  diffuser  regions,  where 
the  second  term  in  the  RHS  of  Eq.  ( 1 6)  is  zero.  Similarly  there 
is  no  solid-state  current  in  the  membrane,  since  the  mem¬ 
brane  is  impermeable  to  electron  flow.  Both  phase  currents 
are  present  in  each  catalyst  layer. 

2.4.  Porous  media  correlations 


The  following  correlations  are  used  to  correct  the  plain 
media  properties  for  porous  media.  The  Rumpf  and  Gupte 
correlation  [  19,20]  is  used  for  permeability.  This  pertains  to 
a  packed  bed  of  spherical  particles  of  average  diameter,  d. 


kp  — 


s5-5d2 

5.6 


(17) 


The  geometric  average  of  the  solid  and  fluid  conductivity  is 
used  to  determine  the  effective  thermal  conductivity  [  19,21], 


ke 

kf 


(18) 


The  dispersion  diffusivity  is  given  by  the  following  relation¬ 
ship,  valid  for  low  Peclet  number  flows  [19,22]: 


Dd 

Off 


=  3(1  — e) 


(19) 


The  total  effective  thermal  diffusivity  is  given  as  [19], 


Dj  _  ke  sDd 
af  kf  af 


(20) 


This  value  depends  on  solid  and  fluid  properties  and  compo¬ 
sitions,  as  well  as  porous  media  characteristics. 

The  effective  gas  diffusivity  is  related  to  the  plain  media 
diffusivity  as  follows  [19,23], 


_  2e 

Dmf  3  -  e 


(21) 


The  electrical  and  ionic  conductivities  are  corrected  for 
porous  media  using  the  Bruggemann  correlation  [24],  where 
fiphase  is  the  volume  fraction  occupied  by  the  phase  through 
which  the  respective  current  flows. 


aeff 


=  as 


1.5 

phase 


(22) 


The  gas  mixture  density  is  found  from  the  ideal  gas  law: 


P  /  nij  \ 

Ur  Ui ) 


(23) 


All  other  fluid  properties,  heat  capacity,  kinematic  viscos¬ 
ity,  thermal  conductivity  and  thermal  diffusivity  are  mass- 
averaged  for  the  gas  mixture.  The  plain  media  gas  pair 


The  unified  (single  domain)  approach  is  used  in  this  model. 
Thus,  boundary  conditions  are  required  at  the  ends  of  the 
domain,  i.e.  at  the  anode  gas  channel/gas  diffuser  interface 
(x  =  0,  boundary  1 )  and  the  cathode  gas  channel/gas  diffuser 
interface  (v  =  0.7  mm,  boundary  2).  All  state  properties  and 
fluxes  are  continuous  at  all  interfaces. 

In  the  experimental  set  up  of  Wang  et  al.  [18],  the  cell 
is  run  at  150  °C  and  1  atm,  and  supply  gases  are  humidi¬ 
fied  at  28  °C.  Thus,  the  total  gas  pressure  at  boundaries  1 
and  2  are  equal  to  the  respective  supply  gas  pressures,  i.e. 
Po  =  101 ,300  Pa.  The  temperature  at  both  boundaries  is  equal 
to  the  respective  supply  gas  temperature,  i.e.  Tq  =  423  K.  The 
partial  pressure  of  water  vapor  is  equal  to  the  vapor  pres¬ 
sure  at  28  °C,  which  is  0.03782  bar.  Thus,  the  mass  fraction 
of  hydrogen  at  boundary  1  is  0.741,  while  the  mass  fraction 
of  oxygen  at  boundary  2  is  0.979.  This  is  the  mass  fraction 
of  oxygen  when  pure  oxygen  is  used  as  oxidant  rather  than 
air.  Since  there  is  no  flux  of  hydrogen  and  oxygen  across  the 
membrane,  the  concentration  gradients  of  each  are  set  to  zero 
at  the  interfaces  of  the  membrane  and  the  respective  catalyst 
layers. 

The  solid  phase  potential  is  arbitrarily  set  to  zero  as 
a  reference  at  boundary  1 ,  and  is  set  to  Ecen  —  Elev  at 
boundary  2,  where  Pcc\\  is  the  desired  cell  potential  and 
£Vev  =  1.154  V  is  the  reversible  cell  potential  at  423  K  and 
1  atm.  The  solid  potential  gradient  is  set  to  zero  are  both  cat¬ 
alyst  layer/membrane  interfaces,  since  there  exists  no  solid 
phase  current  in  the  membrane.  The  electrolyte  phase  poten¬ 
tial  only  exists  in  the  catalyst  layers  and  the  membrane.  Since 
there  is  no  electrolyte  current  in  the  diffuser  regions,  the 
electrolyte  potential  gradient  is  set  to  zero  are  both  catalyst 
layer/diffuser  region  interfaces. 

2.6.  Catalyst  layer  calculations 

In  Wang  et  al.  [18],  E-TEK  electrodes  with  a  Pt  loading 
of  0.5  mg  cm-2  were  used  at  the  anode.  For  the  cathode,  a 
homemade  electrode  was  used,  where  Pt  was  sputtered  on 
the  commercial  electrodes  to  enhance  its  activity.  The  total 
catalyst  loading  at  the  cathode  was  2  mg  cm'2.  Although  not 
reported,  the  catalyst  region  usually  consists  of  30  wt.%  Pt/C 
(f=  0.3).  Since  this  paper  is  interested  in  the  effect  of  porous 
media  characteristics  on  fuel  cell  performance,  it  is  necessary 
to  perform  catalyst  region  calculations. 

The  catalyst  region  is  a  fused  layer  between  the  membrane 
and  the  electrode  sputtered  with  catalyst  particles.  In  this 
region,  there  exists  membrane  phase,  diffuser  solid  and  void 
regions,  as  well  as  catalyst  particles.  There  must  be  adequate 
physical  contact  among  all  the  phases  for  effective  electro¬ 
chemistry  (see  Fig.  1).  The  catalyst  region  consists  of  the 
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following  regions  in  the  respective  proportions: 

membrane  phase:  em  c; 
diffuser  solid  region:  (1  —  em,c)  (1  —  e®); 
diffuser  void  regional  —  emjC)  (e(d))  —  mQ/p0  Zc; 
Pt/C  catalyst  particles:  mjp0  lc. 


where  m0  and  pQ  are  the  average  catalyst  particle  mass  load¬ 
ing  per  MEA  cross-section  and  density,  respectively.  If  /«pt 
is  the  mass  loading  of  Pt  per  MEA  cross-section,  and/is  the 
fraction  of  Pt  in  the  Pt/C  particles,  then 


m0 


Po  =  — 


m  Pt 

/ 

f_ 

PPt 


1  -/ 
Pc 


-l 


(24) 

(25) 


Thus,  the  porosity  of  the  catalyst  region  can  be  determined  in 
terms  of  the  porosity  of  gas  diffusion  region  and  the  catalyst 
packing  properties. 


e(  o  =  £(d)(1_emc)_^ 

*C  J 


1  -/ 
Pc 


(26) 


It  can  also  be  shown  that  the  effective  Pt  surface  area, 


6m  pt 
k  PPt  <7pt 


(27) 


The  diameter  of  the  Pt  particles  is  typically  20  A  [25].  So,  the 
effective  surface  area  can  be  increased  either  by  increasing 
the  catalyst  mass  loading  or  decreasing  the  mean  catalyst 
particle  size. 

Because  the  Pt  loading  is  different  at  the  anode  than  at  the 
cathode,  the/value  would  also  be  different.  It  can  be  shown 
that  at  the  cathode,  the/value  is  given  as, 


fc  = 


m  pt,c 


m  Pt,  e  +  «*Pt, a 


(28) 


3.  Results  and  discussion 

The  governing  equations  are  solved  using  FEMLAB®  3.0. 

The  domain  is  divided  into  768  finite  elements — 288  in  each 

gas  diffuser  sub-domain,  48  in  each  catalyst  layer,  and  96  in 
the  membrane.  A  finite-element  method  is  used  to  solve  the 
system  of  coupled  partial  differential  equations. 

Tables  2-4  give  a  list  of  numerical  values  used  in  the 

computations.  This  will  be  called  the  base  case — the  same 

operating  and  geometric  conditions  as  the  experimental  set 

up.  From  the  paper  by  Wang  et  al.  [18],  the  ionic  conductivity 
of  the  80  p.m  thick  membrane  is  1.87  S  nr  1 . 

However,  they  did  not  report  the  exchange  current  densi¬ 

ties  used  in  their  work.  For  MEAs  with  a  Nafion®  membrane, 

much  investigation  has  been  conducted  into  the  factors  affect¬ 

ing  the  electrode  activity,  thus,  the  exchange  current  density 

is  well-established.  Such  investigation  has  not  yet  been  con¬ 
ducted  for  the  PBI  membrane/Pt  catalyst  interface.  This  value 


Table  2 

Membrane  and  diffuser  properties 
Membrane 


Thermal  conductivity,  km  [26]  (W  m-1  K“ 

-1) 

40 

Ionic  conductivity,  am  [18]  (Sm-1) 

1.87 

Thickness,  /m  [18]  (m) 

8  x  10”5 

Diffuser  (graphite) 

Thermal  conductivity,  (W  m-1  K-1) 

1.15 

Electrical  conductivity,  (S  m-  1 ) 

120 

Porosity,  e 

0.4 

Particle  diameter,  d  (m) 

1.2  x  10~4 

Permeability,  kp  (m2) 

1.8  x  10~u 

Thickness,  1&  (m) 

2.6  x  10~4 

Table  3 

Catalyst  layer  properties 

Anode 

Cathode 

Reference  exchange  current  density,  ij?f 

6 

4  x  10-10 * * * 

(Am-2) 

Effective  catalyst  surface  area,  a  (m-1) 

1.4  x  107 

5.6  x  107 

Transfer  coefficient,  a 

0.5 

2 

Concentration  parameter,  y 

0.5 

1 

Pt  mass  loading,  mpt  [18]  (kgm-2) 

0.005 

0.02 

Weight  fraction  of  Pt  in  Pt/C,/ 

0.3 

0.63 

Porosity,  s 

0.106 

0.092 

Permeability,  kp  (m2) 

1.2  x  10"14 

5.5  x  10~15 

Reference  pressure,  Pref  (Pa) 

1.013  x  10s 

1.013  x  105 

Pt  density,  ppt  (kgm-3) 

21400 

21400 

C  density,  pc  (kgm-3) 

1800 

1800 

Fraction  of  membrane  phase  in  catalyst 

0.4 

0.4 

layer,  £m,c 

Catalyst  layer  thickness,  lc  (m) 

5  x  10~5 

5  x  10“5 

Pt  particles  diameter,  r/pt  (m) 

2  x  lO^9 

2  x  10-9 

depends  on  a  number  of  factors  including  the  manufacturing 
and  assembling  techniques  used  to  form  the  MEA.  Since  the 
exchange  current  density  has  not  been  reported,  in  our  model, 
we  adjust  the  cathode  exchange  current  density  to  best  fit  the 

experimental  data. 

It  is  found  that 

an  exchange  current 

density  of 

4  x  10  10  Am  2  at  the  cathode  and  6  Am  2 

at  the  anode 

give  the  best  fit.  Thus,  these  values  will  be  used  as  the  ref¬ 
erence  exchange  current  densities  at  423  K  and  1  atm.  Fig.  2 
shows  the  fit  between  the  model  predictions  and  the  experi- 

mental  polarization  data. 

Table  4 

Fluid  properties  [27,28] 

Oxygen 

Water  vapor 

Hydrogen 

Molar  mass,  M  (kg  mol- 1 ) 

32  x  10~3 

18  x  10“3 

2  x  10~3 

Thermal  conductivity,  k 

0.0363 

0.03 

0.239 

(Wn-r1  KT1) 

Heat  capacity,  c 

956 

1980 

14450 

(Jkg-1  K-1) 

Kinematic  viscosity,  v 

32.9  x  10“6 

31.11  x  10"6 

302  x  10~6 

(m2  s"1) 

Thermal  diffusivity,  a 

44.4  x  10“6 

30.8  x  10~6 

217  x  10~6 

(m2  s"1) 

Gas  pair  diffusivity,  Dm f 

42.0  x  10“6 

144.0  x  10"6 

(m2  s"1) 
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Current  Density,  I  (A  m'2) 

Fig.  2.  Experimental  and  model  polarization  data. 

Fig.  3  shows  the  potential  variation  within  the  cell  at  lim¬ 
iting  conditions  for  the  base  case.  It  shows  a  0.372  V  ohmic 
overpotential  across  the  membrane  and  a  0.485  V  activation 
overpotential  ( [<pe  —  <ps)  at  the  cathode.  The  former  is  due  to 
low  membrane  conductivity  while  the  latter  is  due  to  slow 
electrode  kinetics  at  the  cathode.  The  ohmic  potential  drop 
across  the  diffuser  region  is  much  lower  (0.090  V  total)  due 
to  its  higher  conductivity  than  the  membrane.  There  is  also  a 
0. 194  V  ohmic  potential  drop  across  the  catalyst  region  since 
this  region  has  a  very  low  effective  ionic  conductivity. 

To  observe  the  effect  of  porous  media  characteristics  on 
the  fuel  cell  performance,  four  factors  are  varied — the  gas 
diffuser  porosity  (ed),  the  anode  /  value  (fa),  a;o,ref  at  the 
cathode,  and  the  membrane  conductivity  (erm).  Fig.  4  shows 
the  effect  of  varying  the  gas  diffuser  porosity  from  its  base 
value.  Although  the  variation  in  performance  is  very  small, 
the  fuel  cell  performs  better  for  a  gas  diffuser  porosity  of 
0.4  than  for  0.3  and  0.5.  This  suggests  that  for  a  given  set 
of  operating  and  geometric  conditions,  there  exists  an  opti¬ 
mum  gas  diffuser  porosity.  The  reason  for  this  is  a  duality 


Current  Density,  I  (A  m'2) 


Fig.  4.  Effect  of  gas  diffuser  porosity  on  polarization  performance. 

between  effective  conductivity  and  permeability  as  the  poros¬ 
ity  changes.  Higher  porosity  implies  higher  permeability  and 
thus,  better  flow  of  reacting  gases.  However,  higher  poros¬ 
ity  means  less  solid  region,  hence,  lower  effective  thermal 
and  electrical  conductivities.  Thus,  there  will  be  larger  ohmic 
potential  drops,  greater  heat  generation  and  higher  tempera¬ 
ture  gradients.  So,  there  is  an  optimum  porosity,  which  gives 
the  best  balance  between  the  two. 

The  IV  curve  for  e  =  0.3  is  the  only  one,  where  a  dis¬ 
tinct  concentration  overpotential  region  is  seen.  Since  there 
is  no  liquid  water,  there  is  no  blocking  of  the  gas  pores  hence 
reduced  transport  limitation.  However,  the  catalyst  poros¬ 
ity  depends  on  the  gas  diffuser  porosity,  and  as  porosity 
decreases,  it  becomes  increasingly  difficult  for  the  gases  to 
permeate  the  catalyst  regions.  When  the  gas  diffuser  poros¬ 
ity  drops  to  0.3  with  the  catalyst  loading  unaffected,  the 
anode  and  cathode  catalyst  layer  porosities  drop  to  0.05  and 
0.03,  respectively  (Eq.  (26)).  Thus,  it  is  not  surprising  that 
mass  limitation  becomes  increasingly  significant  as  porosity 
decreases. 
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Fig.  3.  Overpotentials  at  limiting  conditions  for  the  base  case. 
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Fig.  5.  Effect  of  Pt  weight  fraction  of  polarization  performance. 

Fig.  5  shows  the  effect  of  varying  fd.  Increasing /means 
that  the  carbon  content  in  Pt/C  is  reduced  while  the  Pt  loading 
remains  unaltered.  This  results  in  a  less  packed  (more  porous) 
but  equally  active  catalyst  region.  Note  that  according  to  Eq. 
(28),  varying /a  has  a  corresponding  effect  on/c.  It  must  be 
noted  though  that  a  certain  amount  of  carbon  is  required  to 
house  the  Pt  particles;  so, /cannot  approach  unity.  Neverthe¬ 
less  the  IV  curves  are  barely  distinguishable  for/a  varying 
from  0.25  to  0.6.  However,  for/a  =  0.22,  there  is  a  sharp  drop 
in  IV  performance.  The  reason  for  this  is  that  for  the  base 
case  condition  (f-d  =  0.30),  the  catalyst  regions  are  not  heavily 
packed  (s  ~  0. 1 ).  As/a  decreases,  the  catalyst  layer  poros¬ 
ity  becomes  lower  hence  mass  transport  limitations  become 
more  significant. 

For/a  =  0.22,  the  anode  and  cathode  catalyst  layer  porosi¬ 
ties  are  0.04  and  0.02,  respectively  (Eq.  (26)).  This  would 
explain  the  concentration  overpotential  region  in  the  IV  curve 
for  /a  =  0.22.  Fig.  6  shows  the  oxygen  concentration  for 


Fig.  7.  Effect  of  cathode  catalyst  activity  on  polarization  performance. 

/a  =  0.22  and  0.30  at  limiting  conditions  (£Ceii  =0)-  The  fact 
that  the  mass  fraction  quickly  drops  to  zero  for  /,  =  0.22  sug¬ 
gests  that  oxygen  is  not  able  to  permeate  the  entire  catalyst 
region,  resulting  in  mass  limitation  effects  and  implying  inef¬ 
ficient  catalyst  utilization.  As/a  increases,  the  porosity  of  the 
catalyst  layer  increases,  hence  it  becomes  more  permeable  to 
the  flow  of  reactant  gases.  As  long  as  the  catalyst  porosity  is 
not  very  close  to  zero,  the  fuel  cell  performs  independently 
Of/ 

Fig.  7  shows  the  IV  curves  for  different  values  of  a;o,ref- 
This  value  can  be  improved  either  by  increasing  the  effec¬ 
tive  Pt  surface  area  or  by  increasing  the  exchange  current 
density.  The  latter  can  be  achieved  by  better  manufacturing 
and  assembling  methods,  which  ensure  better  physical  and 
chemical  contact  between  participating  phases  in  the  catalyst 
region.  Predictably,  the  fuel  cell  performance  increases  sig¬ 
nificantly  as  the  catalyst  activity  increases.  The  limiting  cur¬ 
rent  approximately  doubles  as  a;o,ref  increases  from  0.0045 


x(m) 


Fig.  6.  Oxygen  concentration  across  the  catalyst  layer  at  limiting  conditions. 
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Fig.  8.  Effect  of  membrane  conductivity  on  polarization  performance. 

to  4.5Am^.  Note  that  this  is  not  an  unrealistic  improve¬ 
ment  given  that  the  corresponding  value  for  an  MEA  using  a 
Nation®  membrane  is  500  Am-3  [29]. 

Fig.  8  shows  the  effect  of  increasing  the  membrane  con¬ 
ductivity.  The  IV  curves  show  a  drastic  improvement  as  the 
ionic  conductivity  approaches  that  of  Nation®.  The  value  of 
9.6  S  m_  1  is  the  maximum  ionic  conductivity  reported  for  a 
composite  of  PBI  in  the  literature.  The  ionic  conductivity  of 
Nation®  is  17  S  nr1. 

So,  Figs.  7  and  8  show  that  there  is  tremendous  poten¬ 
tial  for  improvement  in  PBIFC  performance  by  increasing 
the  membrane  conductivity  and  enhancing  the  catalyst  activ¬ 
ity.  Based  on  the  above  results,  an  ideal  performance  can  be 
conceptualized  for  a  PBIFC.  Consider  and  MEA  with  cat¬ 
alyst  particles  containing  60wt.%  Pt/C,  fl;o,ref  =  4.5  Am-3 
(the  highest  value  in  Fig.  7),  am  =  9.6  S  nr 1  (the  highest  PBI 
membrane  conductivity  in  the  literature).  All  other  operat¬ 
ing  conditions  are  the  same.  For  these  conditions,  the  opti¬ 
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Fig.  9.  Polarization  and  power  density  for  the  base  case  and  ideal  case. 

mum  gas  diffuser  porosity  is  0.3.  This  represents  the  ideal 
case. 

Fig.  9  compares  the  base  case  performance  with  the  ideal 
performance.  There  is  a  four-fold  increase  in  limiting  current 
density  and  a  five-fold  increase  in  peak  power  density.  The 
peak  power  densities  for  the  base  case  and  ideal  case  are, 
respectively,  1.50  and  7.66  kW  m_2.The  peak  power  density 
for  the  ideal  case  is  in  the  same  order  as  that  of  a  Nation®  fuel 
cell.  A  major  advantage  with  the  intermediate  temperature 
fuel  cell  is  the  absence  of  any  distinguishable  concentration 
overpotential  region  in  the  IV  curve,  implying  that  a  high 
power  density  is  sustainable  over  a  wide  range  of  current 
densities.  With  no  need  for  water  input  into  the  fuel  cell,  and 
sustained  power  over  a  wide  range  of  loading  conditions,  this 
makes  the  PBIFC  very  attractive  for  automotive  applications. 

Fig.  10  compares  the  temperature  distribution  across  the 
MEA  for  the  ideal  case  and  the  base  case.  For  both  cases,  the 
temperature  profile  is  selected  at  the  point  of  optimum  power 


Fig.  10.  Temperature  distribution  at  optimum  power  density  for  the  base  case  and  ideal  case. 
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(i.e.  at  a  cell  voltage  of  0.4  V  for  both  cases).  The  ideal  case 
shows  a  maximum  temperature  rise  of  3.2  K  compared  with 
0.8  K  for  the  base  case.  The  reason  for  this  difference  is  that 
for  the  ideal  case,  the  current  density  is  much  higher  resulting 
in  greater  ohmic  heating  and  heat  of  reaction.  For  the  ideal 
case  at  limiting  conditions,  the  maximum  temperature  rise  in 
the  cell  is  9  K  (not  shown).  The  typical  temperature  rise  for  a 
Nafion®  fuel  cell  is  2-3  K  [30-32].  However,  heat  removal  is 
not  as  critical  for  PB1  as  it  is  for  Nafion® .  PBI  is  already  more 
stable  at  higher  temperatures,  but  most  importantly,  there  is 
no  danger  of  dehydration  of  the  PBI  membrane  since  it  does 
not  require  liquid  water. 

4.  Conclusions 

A  mathematical  model  of  a  PBIFC  was  developed  with 
predicts  the  polarization  characteristics  as  well  as  the  one¬ 
dimensional  transport  phenomena  across  the  MEA.  Model 
results  compare  well  with  experimental  data.  Parametric  stud¬ 
ies  into  the  porous  media  characteristics  of  the  fuel  cell  show 
that  the  fuel  cell  performance  is  good  as  long  as  the  catalyst 
regions  are  sufficiently  porous.  Mass  limitation  effects  are 
only  present  when  the  catalyst  regions  are  heavily  packed, 
resulting  in  extremely  low  permeability  to  reactant  flow. 

The  model  results  also  predict  that  substantial  improve¬ 
ments  in  fuel  cell  performance  are  possible  as  the  conduc¬ 
tivity  of  the  PBI  membrane  and  the  cathode  catalyst  activity 
increase.  Research  is  needed  to  better  understand  the  inter¬ 
facing  between  the  PBI  membrane  and  the  Pt  catalyst,  and 
to  find  ways  to  make  the  catalyst  more  active.  There  exists 
great  scope  for  research  into  finding  composites  of  PBI  whose 
conductivity  approach  that  of  Nafion®,  as  well  as  developing 
techniques  to  better  apply  the  catalyst  loading  and  assemble 
the  MEA.  The  right  combination  of  the  above  may  prove  to 
be  a  significant  breakthrough  in  making  fuel  cell  technology 
commercially  viable. 
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